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ABSTRACT 

Motivated by recent developments in inflationary cosmology indicating the possibility 
of obtaining genuinely open universes in some models, we compare the predictions 
of cold dark matter (CDM) models in open universes with a variety of observational 
information. The spectrum of the primordial curvature perturbation is taken to be 
scale invariant (spectral index n — 1), corresponding to a flat inflationary potential. 
We allow arbitrary variation of the density parameter flo and the Hubble parameter 
h, and take full account of the baryon content assuming standard nucleosynthesis. We 
normalize the power spectrum using the recent analysis of the two year COBE DMR 
data by Gorski et al. We then consider a variety of observations, namely the galaxy 
correlation function, bulk flows, the abundance of galaxy clusters and the abundance 
of damped Lyman alpha systems. For the last two of these, we provide a new treatment 
appropriate to open universes. We find that, if one allows an arbitrary h, then a good 
fit is available for any f2o greater than 0.35, though for Slo close to 1 the required h is 
alarmingly low. Models with Slo < 0.35 seem unable to fit observations while keeping 
the universe over 10 Gyr old; this limit is somewhat higher than that appearing in 
the literature thus far. If one assumes a value of h > 0.6, as favoured by recent 
measurements, concordance with the data is only possible for the narrow range 0.35 < 
flo < 0.55. We have also investigated n ^ 1; the extra freedom naturally widens the 
allowed parameter region. Assuming a range 0.9 < n < 1.1, the allowed range of CIq 
assuming h > 0.6 is at most 0.30 < f2o < 0.60. 

Key words: cosmology: theory - dark matter. 



1 INTRODUCTION 

Even before the announcement of the detection of mi- 
crowave background an isotropics by the D MR experiment 
on the COBE satellite (Smoot et al. 1992), it was realized 
that structure formation models based on cold dark mat- 
ter (CDM) and a flat spectrum of primordial perturbations 
fared considerably better against the data if the matter den- 
sity was reduced by a factor of around three. Most stud- 
ies of this possibility i nvoked a cosmological constant to 
restor e spatial flatness ([Efstathiou, Sutherland fc Maddox 



199C|; [Kofman, Gnedin fc Bahcall 1993Q , with little atten- 
tion being directed to the possibility that the cosmological 
constant may be redundant and the low-density model im- 
plemented in a genuinely open universe. This produces the 
same shape of perturbation spectrum on scales well below 
the curvature radius, but a different normalization and red- 
shift dependence. 

The reluctance to study such models (though general 
arguments in favo ur of an open unive rse w ere developed, 
e.g. Coles & Ellis (1994) and Primack (1995)) arose from a 
widespread belief that inflation, the most plausible candi- 



date for generating the initial density perturbations, could 
give rise to an open universe only in exceptionally fine-tuned 
circumstances. However, open universe inflation models have 
received renewed interest recently, and in particular atten- 
tion has been drawn (Sasaki et al. 1993a, 1993b; Tanaka & 
Sasaki 1994; Bucher, Goldhaber & Turok 1995; Yamamoto, 
Sasaki & Tanaka 1995a; Sasaki, Tanaka & Yamamoto 1995; 
Yamamoto, Tanaka & Sasaki 19 95b; Bucher fc Turok 1995) 



to 



Gott 1982 



the bubb le nucleation model dColcman fc de Luccia 198C 



1 



uth fc Weinberg 1983|: 



Gott & Statler 1984 



Linde 1995; |Amendola, Baccigalupi fc Occhionero 1995). In 



contrast with the situat i on for ordinary model s of inflation 
( [Lyth fc Stewart 1990a} |Ratra fc Peebles 199E| ), this model 
predicts the present value of the density parameter in terms 
of the scalar field potential, without any reference to initial 
conditions. The price that one pays for this is a non-generic 
scalar field potential, which will be even more difficult than 
usual to realize in the context of a sensible particle physics 
model. 

In order to compare structure formation models based 
on open universe inflation models with observational data, it 
is crucial to be able to normalize the amplitude of the power 
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spectrum to the COBE observa tions of cosmic mi crowave 
background (CMB) fluctuations (Bennett et al. 1994), which 



are by far the most accurate available. Anisotropy calcula- 
tions in an open universe present many technical difficul- 
ties and progress to the result has consequently been slow. 
However, an accurate normalization is now available through 
the work of Gorski et al. (1995, henceforth GRSB ), super- 
seding earlier version s by Ratra & Peebles (1994) and by 
Kamionkowski et al. (1994). 



The viability of the open C DM m odel has recently been 
investigated by Ratra & Peebles ( 1994 ) . This paper predated 
the improved normalizations of inflationary open models to 
COBE supplied first by Kamionkowski et al. (1994) and then 
more accurately again by GRSB (see also White & Bunn 
1995). Each of those collaborations provided only a brief ac- 
count of the model's status against observations other than 
those of the microwave background. It is our aim in this pa- 
per to make a more extensive comparison of the model with 
observations. 



2 THE OPEN UNIVERSE POWER 
SPECTRUM 

The model is defined by giving the spectrum V(k) of the 
density contrast 8. Our approach towards constraining the 
model is to utilize linear perturbation theory, applied across 
as wide a range of scales as possible. By considering the for- 
mation of objects such as quasars and damped Lyman alpha 
systems at moderate redshift, it is possible to impose con- 
straints on the spectrum at scales down to one megaparsec 
or less, while COBE probes scales of several thousands of 
megaparsecs, up to and even above the curvature scale. Be- 
tween these extremes, a variety of different constraints can 
be applied. 

All of the observations except the CMB anisotropy 
probe scales which are small compared with the Hubble dis- 
tance, so we can use the Newtonian description of density 
perturbations to describe them. At any epoch well after mat- 
ter do mination oota in, tho power opoctrum of tho donoity 
contrasr is 



V(k) 



(i) 



Here a is the scale factor of the universe, H — a/ a is the 
Hubble parameter (dots signifying time derivatives), k is the 
comoving wavenumber (in inverse megaparsecs) and we are 
defining V{k) as the power per unit logarithmic interval of 
k. The transfer function T(k) specifies the scale-dependent 
effect of the evolution of the density perturbation between 
horizon entry and matter domination, and is normalized to 
unity on large scales. The factor g(0.) is introduced to allow 
for the growth law for perturbations in an open universe. It 
gives the total suppression of growth in an open universe rel- 
ative to a critical densit y universe, and is accurately g iven^ 
by the fitting function (Carroll, Press fe Turner 1992|) 



g(Q) = -a 



(2) 



Numerical tests indicate that this fitting function is accurate 
to within one per cent for Hq of interest. 



Finally, the quantity Sh specifies the overall normalization 
of the present-day spectrum. Its independence of k indicates 
the assumption of an inflationary model leading to a scale- 
invariant spectrum; in typical infl ationary models one w ould 
expect some deviation from this (Liddle & Lyth 1993) and 



we shall consider this at the end of the paper. Inflationary 
models also predict the existence of gravitational waves, but 
their effect on the COBE normalization has not yet been 
successfully calculated and so we assume they are negligible. 
The value of 5h, when fixed by the COBE observations as 
discussed below, depends on Qq, but it has only an extremely 
weak dependence on Ho which can be comfortably ignored 
(throughout, a subscript denotes the present day). 

Many observations do not allow one to impose con- 
straints on the power spectrum itself, but instead place 
constraints on the dispersion of the density contrast a(R) 
smoothed on a comoving scale R. We shall always use a 
top-hat filter W(kR) defined by 



W{kR) 



f sin(kR) cos(kR) \ 
" " ~JkRj 3 {kRf ) 



(3) 



to perform the smoothing. The dispersion of the smoothed 
density contrast is easily calculated from a theoretical power 
spectrum as 



\R) = 



Ah 

V(k) W 2 {kR) ^ . 

k 



(4) 



The prediction for the abundance of objects of various types 
has a very simple interpretation as a constraint on u(R), 
which cannot be easily represented as a power spectrum con- 
straint. 

Before proceeding to a full account of the data and their 
interpretation, let us be more specific regarding our assump- 
tions. The parameters that we shall consider as freely vari- 
able are the total present density fio and the present Hub- 
ble parameter h (in units of 100 kms _1 Mpc _1 ). An im- 
portant contribution to be taken into account is the bary- 
onic component of the density, £Ib, which we take to be 



fixed by nucleosynthesis as Qsh — 0.0125 (Walker et al 



199l| )t. In the presence of baryons, the usual scaling law of 



the transfer function with Qoh (which is exact only for zero 
baryon density) can be replaced by an empirical scaling law 
with Qph e xp(— S Ib — Qb/^o)- This law was discovered by 
Sugiyama (1994), and gene ralizes a scaling law advertised 
by Peacock & Dodds (1994) to the case where Qo < 1- Al- 
though Sugiyama's calculations were made for the case of 
a flat universe with a cosmological constant, the difference 
between that and the present case only sets in long after 
the universe is matter dominated and so the shape is the 
same in our case. The different overall normalization of the 
spectrum between the two cases is of course included in the 
COBE normalization that we shall carry out. 



We use the transfer function from Bardeen et al. (1986) 



7cDM(g) = 



m(l + 2.34g) 
2.34(7 



(5) 



t We no te th at the more recent analysis of Copi, Schramm & 
Turner (1995) suggests that the traditional upper limit from 



Walker et al. (1991) may be relaxed somewhat, though not suffi- 
ciently to impact on our results. 
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41 -1/4 



[l + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71g) 4 ] 

with q = k/hF, where the so-called 'shape parameter' T is 
defined as 



r = fio/iexp(— f2a — fiFj/fio) 



(6) 



in accordance with Sugiyama (L994) as discussed aboveQ 

Although h is in principle freely variable, it is deter- 
mined at some level of accuracy by the requirement of a 
reasonable fit to the galaxy correlation function (see further 
discussion below), which demands that F should lie in the 
range [0.22, 0.29] at 95 per cent c onfidence level assumin g 
a scale-invariant power spectrum (Peacock & Dodds 1994). 



Note though that, as fio tends to 1, the required value of h 
to achieve this begins to get uncomfortably small. 

Particularly for high h, one is in danger of a conflict 
between measured ages of stellar populations and the age of 
the universe. In an open universe the age is given by 



to = 



1 



fio -Ho 



fio 



05 



i - n 2 (i - n ( 



\3/2 



cosh 



.! / 2-fi \ 
V fio J 



■(7) 



If one fixes Ho, then higher ages are achieved by lower- 
ing fio- However, we have written it this way to emphasize 
an alternative view, which is that the galaxy correlation 
function more or less fixes (ignoring for now the baryonic 
corrections) the combination fio-ffo- Then the quantity in 
square brackets in the above formula is actually an increas- 
ing function of fio, peaking at 2/3 when fio = 1. Conse- 
quently, at fixed F, the desire for a large age favours a larger 
value of fio- To make this concrete, then fixing F = 0.25 
taking the baryons into account gives the sample values 
fio = 0.2 => h = 1.3 => Age = 6 Gyr; Q = 0.3 => h = 
0.89 => Age = 9 Gyr; fi = 0.4 => h = 0.69 Age = 11 
Gyr; fi = 1.0 =>• h = 0.32 =>■ Age = 20 Gyr. In each case 
the 15 per cent or so uncertainty in F contributes a similar 
uncertainty to the age. 

We shall adopt the extremely conservative view that 
the age should exceed 10 Gyr, though there are many indica- 
tions that the Universe is older (e.g. Demarque, Deliyannis & 
Sarajedini 1991; Stockton, Kellogg & Ridgway 1995) which 
one can use to constrain cosmological parameters without 
reference to large scale structure. 



3 NORMALIZATION TO CO BE 

The most crucial piece of data is the overall normaliza- 
tion of the density perturbation spectra, which we choose 
to match the microwave anisotropics at large angular scales 
measured by the DMR experiment on t he COBE satellite 



(Bennett et al. 1994; Gorski et al. 1994). In the language 
of the usual spherical harmonic decomposition, COBE mea- 
sures the multipoles with I 30, and for a given fio the 
distance subtended at the surface of last scattering is com- 
parable with the curvature if I £ 2yl — fio /fio- With the 
possible exception of the super-curvature modes defined be- 
low, this criterion gives an upper bound on the range of I 



for which curvature can be significant. It may however be 
a considerable overestimate, because for fio < 1 the dom- 
inant contribution to the CMB anisotropy can come from 
distances far closer than the surface of last scattering. In 
any case it allows curvature to affect only I ^ 6 even for 
fio as low as 0.3, which means that at most the lowest few 
multipoles of COBE are likely to be sensitive to curvature. 

To investigate the effect of curvature quantitatively, one 
must first ask how the Newtonian expression (^) should be 
continued t o large r scales. As discussed in detail in Lyth & 
Woszczyna ( 1995 ), a number of issues have to be addressed. 

First, in order to define the density contrast one has to 
specify a slicing of space-time into spatial hyp ersurf aces. We 
make the usual choice that the hypersurfaces are orthogonal 
to comoving observers, corresponding to what is called the 
'gauge invariant' density perturbation. In the era well after 
matter domination (which is the only one that concerns us) 
this is t he same as the 'synchr onous gauge' density pertur- 
bation (|Lyth fe Stewart 1990rj). 



Secondly there is the definition of the spectrum. In dis- 
cussing the stochastic properties of a given perturbation /, 
one assumes that it is a typical realization of a random field 
(an ensemble of functions together with a probability distri- 
bution for them). In both flat and curved space, the spec- 
trum is defined with reference to an expansion in terms of 
eigenfunctions of the Laplacian, being the ensemble average 
of the m odulus squared of the coefficient. Following Lyth & 
Stewart (1990a), we denote the eigenvalue of the Laplacian 
by — (fc/a) , and normalize the spectrum Vf(k) of a generic 
perturbation / so that it gives the power per unit logarith- 
mic interval of k. (By 'power' we mean the ensemble mean 
square contribution to f 2 , which is independent of position.) 
We are taking the random field to be Gaussian, which means 
that each coefficient has an independent Gaussian probabil- 
ity distribution, whose variance is essentially defined by the 
spectrum. (To make this statement precise one needs to take 
account of the fact that k is a continuous, not a discrete, 
variable.) 

Thirdly, there is the range of k over which the spectrum 
is non-zero. If fc _1 is measured in units of the curvature scale 
Hq 1 j\J\ — fio, then it is known that the most general square 
integrable junction can be constructed using only the eigen- 
functions with k 2 > 1. For this reason, cosmologists have 
always assumed that the same is true for the most general 
Gaussian random field. That is, they have assumed that such 
a field can always be generated by keeping only the the sub- 
curvature modes (those with < k~ 2 < 1) as distinct from 
the super-curvature mod es (those with k~ 2 > 1). It has re 



cently been pointed out (Lyth & Woszczyna 1995) that this 
is not so; rather, mathematicians have known for half a cen- 
tury that in order to construct the most general Gaussian 
random field the spectrum (and therefore the eigenfunction 
expansion) needs to run over the full range k 2 > 0. Such 
super-curvature mode s can arise in the single-b ubble mod- 
els of open inflation (Yamamoto et al. 1995b), though as 
we shall see it is a reasonable working hypothesis to assume 
that the effect of these is negligible.^] 



t Note that Peacock &: Dodds ( 1994 ) have a typographical error 
in the transfer function; we have confirmed that their results apply 
to the correct form. 



S A smooth continuation of the spectrum into the super-curvature 
regime would not have a sig nificant effect on the CMB anisotropy 
(Lyth fc Woszczyna 1995j), One can also show that a delta 
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In most of the cosmology literature a different normal- 
ization of the spectrum is adopted, which is denoted by Pf 
rather than by Pf. In flat space, Pf is normalized so that 
k 3 Pf/(2n 2 ) is the power per unit logarithmic interval of k. 
Because super-curvature modes were never considered, this 
definition is customarily generalized to make q 3 Pf / (2tv 2 ) the 



is only 2. In the single-bubble case, there is also the possibil - 



power per unit logarithmic interval of q, where 



k 



1. 



(The motivation for considering q 2 instead of k 2 is that its 
range is q 2 



2n 2 



Vf(k(q))- 



M. 



These preliminaries having been addressed, we are 
ready t o ask what is the correct continuation to large scales 



of the flat space expression (1)? A natural choice is to 



keep equation (|lj) as it stands, either retaining or drop- 
ping the super-curvature modes. If super-curvature modes 
are dropped there are other natural choices, based on the 
alternatively defined spectrum P that we have just dis- 
cussed. One can take P oc q (the usual choice until re- 
cently) or P oc k; these choices multiply e quation (111) by 



(q/k) and q/k respectively. Another choice (Kamionkowski 



Sper ;el 1994 ) , relating to the density contrast smoothed 



over a sphere of variable radius, makes Pag for q < 1 
going smoothly over to P oc q for q <; 1. 

A further possibility is that, instead of focusing on the 
density contrast 5, one can f ocus on the primordial curvature 
perturbation 1Z, given by (Bardeen 1980; Lyth & Stewart 
1990a[)] 



1Z : 



a 2 H 2 
3 + fc 2 



(9) 



On small scales equation (|l]) corresponds to a scale-indep- 
endent spectrum V-n- However, if Viz is taken to be scale 
independent on large scales also, equation (jl|) is multiplied 
by a factor [(3 + k 2 )/k 2 } 2 . At k 2 = 2 this is a factor ~ 5, 
corresponding to a factor Vo~ in the rms perturbation, so 
it is more significant than the ambiguity associated with 
the definition of the spectrum, and the use of q 2 versus k 2 . 
Thus the crucial decision is whether to regard the density 
perturbation or the curvature as the fundamental quantity. 

The usual assumption that the perturbation originates 
as a vacuum fluctuation of the inflaton field decides in favour 
of the curvature, because the inflaton field perturbation S<j> 
is related to t he curvature by 1Z = (H/<j))S<j) which is scale - 



independent ( Lyth fc Stewart 1990a Liddle fc Lyth 1993 ) 
Making the arbitrary assumption of the conformal vacuum 
as the initial state in calculating the inflat on perturbation, 
the s pectrum of 5<f> and theref ore of TZ is flat ( Lyth & Stewart 
1990a ~pi,atra fc Peebles 1995 ). Recently it has been pointed 



out that in the bubble nucleation model the quantum fluctu 
ation of the inflaton field, and hence the spectrum, is calcula- 
ble without recourse to an arbitrary assumption concernin 
the initial vacuum stat e. Acc ording to Bucher et al. (199 
and Bucher & Turok ( 1995 ), Vn varies like coth(7rq). At 
k 2 = 2 (q 2 — 1) this factor is 1.06, and even at q 2 = 0.03 it 



function contribution at k 2 = (the open unive rse Grishchuk- 
Zel'dovich effect) is not compatible with the data ( |Garcfa-Bellido 
et al. 1995D. 



ity of a discrete super-curvature mode (Sasaki et al. 199 



provided the inflaton mass is light enough; however, again 
this should have only a small effect. Hence the difference 
between the conformal vacuum hy pothesis and the bubble 
nucleation scenario is insignificant (Yamamoto et al. 1995a; 
Bucher fc Turok 1995| )fl 

Although present versions of the open inflationary sce- 
nario are not without their problems, the situation does ap- 
pear to have improved recently, in that the fine-tuning of 



initial con ditions required in early models (Lyth & Stew- 



art 1990a) is not necessary in the bubble nucleation model. 



At present it does however still seem necessary to have 
some fine-tuning in the parameters of these models (Linde & 



Mezhlumian 1995). Despite this, open inflation models are 



the natural underpinning of structure formation models in 
an open universe. 

The upshot of the above discussion is that the criterion 
of some smooth continuation suggests a moderate amount of 
ambiguity in the large-scale power spectrum, which is how- 
ever practically eliminated if the perturbation originates as 
a quantum fluctu ation of the inflaton field. According to 
Sugiyama & Silk (1994), even the moderate ambiguity sug- 



gested by smoothness is not very significant, because it af- 
fects only the low CMB multipoles which are poorly deter- 
mined because of cosmic variance. However, in this paper 
we adopt for definiteness the hypothesis of a flat curvature 
spectrum, which is the prediction of inflation. 

Given the spectrum, one can calculate the expected val- 
ues for the multipoles measured by COBE, and compare 
with observations to determine the best-fitting normaliza- 
tion. The calculation is substantially more complex than 
that for critical density models, which is almost analytic, 
and in the literature it has been developed in several stages. 
In this paper we use the most recent and sophisticated de- 
termination, given by GRSB. They take the curvature per- 
turbation spectrum to be flat, and fit the full spectrum 
of anisotropies including the Doppler peak using a method 
based on Fourier analysis on the cut sky for which COBE 
data are available. This normalizati on is m ore sophisticated 
than that of Kamio nkowski et al. (1994| ), who normalized 
to the 10° variance ( Bennett et al. 1994[ ) with a correction 
incorporated for the beam profile and non-orthogonality o f 
the monopole and dipole subtraction (Wright et al. 1994). 
They also arbitrarily increased the error bar to 30 per cent. 

The outcome of the GRSB analysis is that essentially all 
values of tto are capable of providing an acceptable fit to the 
COBE data for a suitable choice of normalization. They do 
not explicitly state the normalization of the power spectrum 
that they get for each £1q ■ However, they do give values 
of cr(8/i _1 Mpc) for specific choices of h, directly calculated 
from their Boltzmann code. We use these to calculate the 
large-scale normalization of the power spectrum (Sh{Qo) in 
equation ([j])), which is independent of h. This can then be 



After this paper was accepted, Yamamoto & Bunn (1995) 



demonstrated explicitly that the normalization of the power spec- 
trum in the single-bubble and conformal vacuum cases is nearly 
identical for any reasonable Qo (though the full details of the fit 
to COBE differ somewhat between the two cases). 
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used to calculate a(R) using equation (^) for any value of h 
by using the appropriate transfer function. 

We find that the normalization can be accurately repre- 
sented, to within 2 per cent for 0.1 < fio < 1, by the fitting 
function 



5fr(n ) = (4.10 + 8.83fi - 8.50fio) x 10" 



(10) 



where we use the GRSB analysis which includes the quadru- 
pole (almost no change arises if the quadrupole is dropped 
from the analysis). 

The normalization from GRSB has an error bar of 8 
per cent (more or less independently of fio), as co mpare d 
with the 30 per cent used by Kamionkowski et al. (1994). 



Although this tighter error bar is certainly more constrain- 
ing, this normalization is quite a bit higher than used by 
Kamionkowski et al. An increase in the normalization gener- 
ally acts in favour of the lower density models when it comes 
to comparing with the observations. As the COBE normal- 
ization has a considerably smaller error bar than other ob- 
servations, on occasion we shall take this normalization as 
fixed, ignoring its error bar. 



4 SMALLER SCALE CONSTRAINTS 

A wide range of observations provide a variety of constraints 
on the power spectrum on scales of order 1 to 100 Mpc. 
These include the distribution of galaxies and clusters, the 
peculiar motions of galaxies and the abundances of various 
objects including clusters, quasars and damped Lyman al- 
pha systems. Our approach is to use only the most powerful 
ones, as described elsewhere for the case fio = 1 (Liddle & 
Lyth 1995; Liddle et al., in preparation). When the spatial 
geometry is changed, all constraints need to be recalculated 
for a variety of reasons, amongst which the primary ones are 
a suppressed rate of perturbation growth at low redshift and 
an amended relation between scale and mass. 



4.1 The galaxy correlation function 

One of the most highly advertised problems with the stan- 
dard cold dark matter scenario is its failure to reproduce 
correctly the shape of the galaxy correlation function on 
scales of tens of megaparsecs, on the reasonable assump- 
tion of a scale-independent bias parameter for galaxies of 
a given type. For CDM models, this is quantified via the 
shape parameter P which we have already introduced, and 
from a detailed analy sis of a compilation of data sets Pea- 
cock & Dodds (1994) obtain the very stringent constraint 
F = 0.255^0 033 at the 95 per cent confidence limit assum- 
ing a scale-invariant power spectrum. Provided one is willing 
to tolerate a sufficiently small h (around 0.35), the shape pa- 
rameter can be fitted in a critical density universe. 

In addition to providing a constraint on the shape pa- 
rameter, the galaxy distribution data also in principle con- 
strain the normalization of the spectrum through redshift 
space distortions and non-linear effects. By choosing a scale 
in the middle of the data the best-fit amplitude can be found 
independently of F; assuming fio = 1 and br = 1, where bi is 
the bias parameter for JiMS-selected galaxies, we find the 
constraint <r(15.1 h~ 1 Mpc) = 0.40 ± 0.03. For general fi , 
Peacock & Dodds provide a best-fitting bias parameter, and 



by fitting for this and processing through the redshift dis- 
tortion factor for general fio one obtains the formal result 
with la error (almost entirely due to the uncertainty in the 
bias) of 



<r(15.1h _1 Mpc) = (0.40 ± 24per cent)/(fi ) , 
where the fitting function /(fio) is given by 
/(fio) = 1.62 + O.Slfio - 2.60fi 2 + 1.31fio . 



(11) 



(12) 



However, the literature contains a widespread range of esti- 
mates of the b ias p arameter (see for example the compila- 
tion in Dekel (1994)), suggesting a true uncertainty larger 
than that advertised by Peacock & Dodds. As this result is 
anyway less constraining than other data, we choose not to 
impose this constraint. 

A chi-squared analysis o f the sixteen data points in ta- 
ble 1 of Peacock & Dodds (1994), taking P and the nor- 



malization as fitting parameters, has 14 degrees of freedom. 
Unfortunately the minimum chi-squared is somewhat low 
(about 12). It is perfectly reasonable that this has occurred 
by chance, though it could also have an origin in weak cor- 
relations of neighbouring data or through non-normal er- 
rors. This prevents us incorporating their full data set into 
a chi-squared analysis along with other data, because such 
an analysis allows other data points to receive a high chi- 
squared in compensation because their data set has so many 
more points. We have tried to evade this by only incorpo- 
rating the shape parameter into a chi-squared test on all the 
data. 



4.2 Bulk flows and POTENT 

For a given present-day amplitude of density perturbations, 
the predicted peculiar velocities depend quite strongly on 
the value of fio, becoming much smaller in the low-density 
case. The best measurements of the bulk flow available are 
those found v ia the POTENT technique of velocity field re- 
constructi on ( Bertschinger fc Dekel 1989 ). For the Mark III 
data set (Dekel 1994), the velocity has been evaluated in 
spheres about our position for a range of radii. However, 
these separate determinations are not independent as the 
rms bulk flow is sensitive to long wavelengths to a much 
greater extent than the density contrast. We therefore con- 
centrate on a single measurement, which is the bulk flow 
smoothed on a scale of 40 hT 1 Mpc. The method used to 
generate this requires an additional Gaussian smoothing on 
12/i _1 Mpc in order to generate the original continuous ve- 
locity field used as a starting point. The theoretical predic- 
tion for the rms bulk flow is therefore given by 



°i{R) = HqQq 2 I W 2 (kR) exp(~ (12 h^k) 2 ) ^^,(13) 



VoAk 
k 2 k 

where W(kR) is the top-hat window given by equation Q 
and the factor fij' 2 is an extremely accurate fitting function 
to the fio-dependent velocity suppression. 

The Mark III POT ENT analysi s gives the bulk flow in 
a 40 ft" 1 Mpc sphere as flDekel 1994]) 



w obs (40 /i _1 Mpc) = 373 ± 50 kms" 



(14) 



and this provides the best estimate of a v (40h x Mpc). The 
error given in expression ([w]) arises from different ways of 



© 0000 RAS, MNRAS 000, 000-000 



6 A. R. Liddle et al. 



dealing with sampling-gradient bias and can thus be thought 
of as reflecting the systematic uncertainty in the POTENT 
analysis. Additionally there is an intrinsic uncertainty in the 
POTENT calculation due to rand om distance errors, which 
at the la level is ~ 15 per cent (Dekel 1994). The obser- 
vational error is dominated by cosmic variance; since the 
mean square bulk flow is the sum of the squares of the three 
velocity components, each of which is Gaussian distributed, 
it follows a x 2 distribution with three degrees of freedom. 
This enables a calculation of the cosmic variance error in 
using the bulk flow as an estimator of the normalization of 
the dispersion of the density contrast, that error being 89 
per cent upwards and 24 per cent downwards at the 68 per 
cent confidence level which notionally corresponds to la. At 
the 95 per cent confidence level the error bars are +273 per 
cent and —43 per cent. We can improve on this by modelling 
the observational errors and convolving with the theoretical 
distribution. Assuming that the error in expression ([lij) cor- 
responds to something like 95 per cent confidence (though 
as it is the smallest error this assumption is insignificant), 
then the convolution of the three types of error results in the 
total error in using the Mark III POTENT bulk flow calcu- 
lation as an estimator of the normalization of the dispersion 
of the density contrast. The increase in the error range as 
compared with cosmic variance alone is not large, the total 
error range being +98 per cent to —25 per cent at the 68 per 
cent confidence level, and +295 per cent to —47 per cent at 
the 95 per cent confidence level. Only the lower limits are 
useful for us. 

We note that a constraint on the value of Jlo can be 
extracted from non-linear effects on the peculiar velocities, 



yiel ding fio <; 0.3 at least at the 2a confidence level (Dekel 
1994) which serves to reinforce our conclusions. 



The scale at which the bulk flow data apply is of order 
1 per cent of the Hubble distance, so one might wonder if 
general relativistic effects might be detectable. The formulae 
that we have given remain valid in that case, provided that 
the density perturbation is defined on hypersurfaces orthog- 
onal to comoving observers, and that the peculiar velo city is 
defined with r espect to worldlines having zero shear (Bruni 
& Lytr ] 1994| ) ■ It is noted in GRSB that, with a different 
choice, the theoretically calculated bulk flow is different by 
several per cent, which is not totally insignificant. This sug- 
gests that a careful analysis of the observations using general 
relativity would be worthwhile, using a specific set of world- 
lines to define the peculiar velocity. In any event, as long 
as Newtonian physics is used to analyse the data there is 
certainly no point in going beyond that framework in the 
theoretical calculation. 



4.3 Object abundances 

In the case of a critical density universe the standard an- 
alytical technique to calculate object abundances relies on 



the use of the Press-Schechter theory (Press & Schechtei 



1974), which has been found thr ough iV-body simu lations 
to provide a good approximation (Lacey & Cole 1994). This 



kind of comparison between analytical techniques and N- 
body simulations has not been performed to the same extent 
for an open universe. However, the derivation of the Press- 
Schechter theory relies solely on statistical arguments; there 
is nothing in it that explicitly relies on the background cos- 



mology. It should therefore also be applicable in an open 
universe. We shall use it to obtain constraints on the abun- 
dances of galaxy clusters and damped Lyman alpha systems. 

Using the Press-Schechter theory, the fraction of the 
matter in the universe that is in collapsed objects above a 
given mass at a redshift z is given simply by 



n(> M(R),z) 



; erfc 



V2a(R, 



(15) 



where S c is the threshold value fixed by comparison with 
JV-body simulations, a(R,z) is the dispersion smoothed on 
scale R at redshift z and 'erfc' is the complementary error 
function. The appropriate value for S c in this expression de- 
pends on the type of collapse one wants to consider, and 
on the type of filter one uses to carry out the smoothing. 
In a critical density universe the spherical collapse of a top- 
hat perturbation is associated with S c — 1.7. Non-spherical 
collapse along all three axes of symmetry is associated with 
higher values for S c , whilst non-spherical collapse along the 
first and second collapsing axes is associated with smaller 
values (e.g. Monaco 1995). As the value of S c is determined 
by the time-scale of collapse of a given type of perturbation, 
one might expect it to be quite sensitive to the background 
cosmology being considered. However, this does not seem to 
be the case when one move s from a critical densit y uni verse 
to an open universe. Lilje (|l992f ), Lacey & Cole (1993) and 
Colafrancesco & Vittorio ( |1994[ ) found that, at least for any 
type of collapse where 5 C < 1.7 in a critical density universe, 
the value of S c varies at most by 5 per cent when one goes 
from an Qo — 1 universe to one with Qo = 0.1. This ap- 
plies at the present epoch, therefore implying that the same 
change in background cosmology will give rise to an even 
smaller variation in S c at higher redshifts, as presently open 
universes approach flatness with increasing redshift. 

The abundance of galaxy clusters is used to constrain 
the present-day power spectrum. In order to constrain 
shorter scales, which are well into the non-linear regime to- 
day, a successful technique is to study objects at high red- 
shift, when those scales were still in the linear regime. By 
using linear theory to scale those constraints to the present, 
one can compare directly with the present-day predicted 
linear power spectrum. The most useful objects on which 
data are available are the damped Lyman alpha systems 
(Lanzctta, Wolfe & Turnshek 1995; Storrie-Lombardi et al. 
1995). These offer a tighter constraint than the quasar abun- 
dance, the latter being weakened by unknown efficiency fac- 
tors such as the required number of generations of quasars, 
and by the uncertainty in the required host galaxy mass. 

We wish to take into account the growth of density 
perturbations between a redshift, say, around four and the 
present. As f2o is decreased, the amount of growth between 
these epochs becomes highly suppressed, which is one of the 
main reasons why the present normalization of the primor- 
dial spectrum is lower than in the critical density case. On 
the other hand, this effect helps with high-redshift object 
formation since, for a given present-day normalization, the 
perturbations at high redshift are substantially higher than 
if the universe were flat. 

In a critical density matter-dominated universe, a(M, z) 
simply grows proportionally to (1 + z) . In an open uni- 
verse, there is a suppression g in growth relative to this 
given by equation (0). This equation can be applied at any 
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epoch, using the redshift dependence of which in a matter- 
dominated universe is given by 



n(z) = n 



1 + n z 



(16) 



One therefore needs to apply the growth factor for a criti- 
cal density universe, correcting for the suppression both at 
the redshift of the observation and at the present, to get a 
constraint on the present-day power spectrum from 



a(M, 2 = 0)= a(M, z) (1 + z) 



gOV) 
s(fi(*)) 



(17) 



4-3.1 Cluster abundance 

A large galaxy cluster has a typical mass of about 10 15 M , 
which corresponds to a linear scale of around 8ft" 1 Mpc. 
Such clusters are relatively rare, indicating that this scale is 
still in the quasi-linear regime. One is then able to use the 
Press-Schechter theory to calculate erg = a(8 fe~ 1 Mpc). To 



our knowledge the first to atte mpt t his was Evrard (1989), 
followed by Henry & Arnaud ( 1991 ), Both these analyses 



were only valid for a critical density CDM universe, and 
though using different observations they reached ess entially 
the same result. Then White, Efstathiou & Frenk ( 1993a ) 
again obtained a result similar to the previous two, and ex- 
tended the analysis to a flat CDM universe with non-zero 
cosmological constant. Our analysis is similar to theirs ex- 
tended to an open universe, the main difference being that 
we attempt to take into account that clusters with equal 
mass which virialize at different redshifts have distinct prop- 
erties, like velocity dispersion and X-ray temperature, at the 
present. 

A variety of different observations are available concern- 
ing the abundance of clusters. To use the Press-Schechter 
theory, it is vital to have good mass estimates as well as 
an estimate of the number density. Galaxy cluster cat- 
alogues assembled through optical selection from photo- 
graphic plates, even disregarding the subjective nature of 
such selection, suffer from possible errors in cluster identifi- 
cation due to foreground and background contamination in 
the galaxy counts. Furthermore, the velocity dispersion, the 
optical observable most directly related to the cluster mass, 
is prone to serious projection effects and possible velocity bi- 
ases. On the other hand, cluster identification through X-ray 
emission is free from foreground and background contamina- 
tion, as X-rays are only produced in deep potential wells, and 
the X-ray observable most directly associated with the clus- 
ter mass, the mean X-ray temperature, is only very weakly 
affected by projection effects. Accordingly, we choose to use 
X-ray instead of optical data. 

The observed number density of clusters per unit tem- 
peratu re, n (kBT), at z = was calculated by Henry & Ar- 



naud (1991). They found that clusters with a mean X-ray 
temperature of 7 keV have a present number density 



n(7 keV, 0) 



2.0±l° x 10" V Mpc" 3 keV" 1 . 



(18) 



The comoving number density of clusters in a mass in- 
terval dM v about virial mass M v at a redshift z is obtained 
by differentiating equation ( |l5| ) with respect to the mass 
and multiplying it by pb /M v , where pb is the comoving back- 
ground density (a constant during matter domination) , thus 



giving 

n(M v ,z)dM v = 



(19) 



Pb 



dA(z) 



7T M v A 2 (z) dM v 



exp 



2A 2 (z) 



dM v , 



where A = cr(ri,) with tl the comoving linear scale asso- 
ciated with M v , r£ = 3M v /47T/9b. Traditionally the cluster 
abundance is used to constrain the dispersion at 8 ft" 1 Mpc, 
and the quantity A is specified by an analytic approximation 
to the power spectrum in the vicinity of this scale. Generally, 
one can write 



8 ft- 1 Mpc/ 

For the CDM spectra we adopt the form (for n = 1) 



7 (r L ) = (0.3r + 0.2) 



2.92 + log 



I'L 



8 ft" 1 Mpc 



(21) 



This is a more sophisticated analytic approximati on tha n 
the power- law approximation used by White et al. (1993a); 
the open universe calculation requires accuracy over a wider 
range of scales (note also that their V has a slightly different 
definition). This approximation is accurate to within 1 per 
cent for n, within a factor of 4 of 8 ft" 1 Mpc for the F values 
in which we are primarily interested. 

Note that, in any CDM model, y(t'l) is redshift inde- 
pendent since the growth of perturbations is independent 
of scale. Using expression (^) to calculate the derivative in 
equation (jl^) we therefore get 



n(M v , z) dM v 



(22) 



2 p b 2.92(0.3F + 0.2)<5 C 



Mv 2 



3A(z) 



■ exp 



2A 2 (z) 



dM v . 



As large clusters are relatively rare, it is reasonable to 
assume that shear did not play an important part during 
their collapse, which to a good approxi mation can then b e 
considered to have occurred spherically (Bernardeau 1994). 
Nevertheless, we shall include an assumed la dispersion of 
±0.1 in the value of 5 C . Bearing in mind that varying the 
background cosmology has a negligible effect on the value of 
5 C we then use <5 C = 1.7 ±0.1 when needed for all our models 
at all z. 

For the type of models we are considering, Hanami 
has shown that 



( [1993D 



M v oc Q" 1/2 x~ 1/2 (l + z c y 3/2 (k B T) 3/2 h' 



where 
A 



l±(fi" ' 8 -l)(l±^- 5 z m )" n o 



(23) 



(24) 



Here z c and z m are the redshifts of cluster virialization and 
turnaround respectively; they are related by the expression 
(1 ± Zm) ~ 2 2 / 3 (l ± z c ). The scalings in equation ( pgj ) have 
been found through hydrodynamical iV-body simulations to 
hold remarkably well in a Qo = 1.0 CDM model (Navarro, 
Frenk & White 1995). 

In order to normalize equation ( p3| ) we use results from 
the hydrodynamical iV-body simulatio ns for a £Iq = 1.0 
CDM model performed by White et al. ( |l993b| ). From a cat- 
alogue of 12 simulated clusters with a wide range of X-ray 
temperatures they estimated that a cluster with a present 
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mean X-ray temperature of 7.5 keV corresponds to a mass 
within one Abell radius (1.5 hT 1 Mpc) of the cluster centre of 
Ma = (1.10±0.22) x 10 15 h' 1 Mq. The error arises from the 
dispersion in the catalogue and is s upposed to represent the 
la significance level. White et al. (1993b) also found that 
the simulated clusters had a density profile in their outer 
regions approximately described by p c {r) oc r~ 2, 4±0,1 . This 
same result wa s obt ained by Metzler & Evrard (1994) and 
Navarro et al. (1995). Bearing in mind that the cluster virial 



radius in a Qq = 1-0 universe encloses a density 178 times 
the background density, it is then straightforward to cal- 
culate the cluster virial mass from Ma. Through a Monte 
Carlo procedure, where we assume the errors in Ma and in 
the exponent of p c (r) to be normally distributed, we find 
M v = (1.23 ± 0.32) x 10 15 hT 1 M for a cluster with a 
present mean X-ray temperature of 7.5 keV in a flo — 1.0 
universe. Assuming that such a cluster virialized at a red- 
shift of z c ~ 0.05±0.05 (e.g Metzler & Evrard 1994; Navarro 
et al. 1995), we can now normalize equation (ps)): 

M v = (1.32 ± 0.34) x 10 15 x (25) 

0-1/2 -1/2 /i i N-3/2 ( k B T \ 3//2 ! 

This result i s in v ery close agreement with the one obtained 
by Evrard ( |l99C| ) from his own hydrodynamical iV-body 
simulations. Hence the virial mass M v for a cluster with 
a present mean X-ray temperature of 7 keV is given by 

(1.2±0.3) x 10 15 fi ( 7 1/2 x~ 1/2 (l+2c)" 3/2 ft" 1 M .(26) 



M v 



As' onc can 3cc from equation (|23 | ), the relation between 

the cluster virial mass and its mean X-ray temperature de- 
pends on the redshift of cluster virialization. One therefore 
expects that at the present there will be some dispersion in 
the virial masses of clusters with the same mean X-ray tem- 
perature. This dispersion increases with decreasing Slo, as, 
due to the slower growth of density perturbations in lower 
Qo models, cluster formation at a given scale proceeds over 
a greater redshift interval. 

According to Press-Schechter theory the comoving 
number density of clusters with virial mass in an interval 
dM v at M v which virialize in a redshift interval dz at red- 



shift z and survive until the present is given by ( sasaki 1994 ) 



N(M V , z) dM v dz = 



(27) 



n(M v ,z) da»(z) 



A 2 (z) a s (z) 



dz 



a s (z) 



a s (z = 0) 



dM v dz , 



where crg(z) and da$(z)/dz are calculated using equation 
(p"7j). In equation ( |27j ) the expression within the square 
brackets gives the formation rate of clusters with virial mass 
M v at redshift z, whereas the fraction outside gives the prob- 
ability of these clusters surviving until the present. If one 
now assumes that at each redshift z the cluster virial mass 
M v in equation ( |27j ) is determined by expression ( p^ ) with 
z c — z, then equation ( p7j ) gives the comoving number den- 
sity of clusters per unit mass which virialize at each red- 
shift z and survive up to the present such that they have a 
mean X-ray temperature of 7 keV at the present. Through 
the chain rule we can then determine the comoving num- 
ber density of clusters per unit temperature that virialize at 
each redshift z and survive up to the present such that they 



have a mean X-ray temperature of 7 keV at the present 

dM v 



N(k B T,z)d(k B T) dz 



d(k B T) 



N(M v ,z)d(k B T) dz, 



= ±—LN(M v ,z)d(k B T)dz,(28) 
2 k B l 

where the second equality uses equation (p3|). We therefore 
have 



N(k B T, z) d(k B T) dz = 
3 M v <5 2 



(29) 



n(M v ,z) dag(z) 



2 k B T A 2 (z) as(z = 0) dz 



d(k B T) dz. 



Numerically integrating this expression from z — to 
z — oo then gives the present comoving number density 
of clusters per unit temperature with a mean X-ray tem- 
perature of 7 keV as a function of Slo and of the present 
value of ag,- Comparing with the observational value given 
by equation ( |l8| ) we thus get <T8(fio). We find to a good 
approximation that 



a 8 = 0.60 



+ 32 per cent 
-24 per cent 



c(n ) 



(30) 



where C(fi ) = 0.37 + 0.13fi - 0.02Slg is a fitting function 
representing the changing power-law index of the Qo depen- 
dence. We have computed the uncertainty using a Monte 
Carlo method; it arises from the dispersions in the obser- 
vational value of r, in the assumed value for S c , and in 
expressions (|l^) and (^6|). The confidence level quoted in 
equation (Bfj) is 95 per cent. Further details of th e calcula- 
tion of the u ncertainty will be provided elsewhere (Viana & 



Liddle 1995) 



4-3.2 Damped Lyman alpha systems 

Many types of model with Qq = 1, such as mixed dark 
matter models, are stro ngly constrained by data on damped 
Lyman alpha systems (jMo fc Miralda-Escude 19 94 , 



p. 



maim fc Chariot 1994; |Ma fc Bertschinger 1994; Klypin et 



Kauff- 



al. 1995). However, the constraint becomes weaker as S7o is 



reduced, as we will now see. 

In stead of the widely quoted data of Lanzetta et al. 
( [l995| ), we use the recent data of Storrie-Lombardi et al. 



( |1995| ) which revises downwards|[J] the estimated abundances 
at a redshift of around 3 and provides a new estimate at 
redshift 4. The strongest constraint comes from the redshift 
4 point, and so we shall concentrate on it. However, the 
constraint is not significantly weakened if the redshift 3 point 
is used instead, and in any case we shall see that these data 
are not very constraining for open CDM models. 

At redshift z — 4, the contribution of the damped Ly- 
man alpha systems to the density in baryons is estimated 



Or 



= (0.0011 ± 0.0002) h~ 



1 + 9,qz 
l + z 



(31) 



where we have conservatively assumed that all the gas in 
these systems is in the neutral state. Remembering that we 



II Note that this still ignores the effect of gravitational lensing, 
which it is claimed can reduce the estimated abundance by a 
further 50 per cent (Bartelmann & Loeb 1995). 
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Figure 1. The constraints plotted in the Qo-h plane, assuming baryon density as given by standard nucleosynthesis. All models are 
normalized to the COBE data as given by GRSB, and the constraints shown are all at 95 per cent confidence. The solid lines are limits 
on the shape parameter. The dashed lines are limits from the cluster abundance. The dot-dashed line is the lower limit from POTENT 
(since it comes from confidence intervals, formally it is 97.5 per cent as a lower limit, the upper limit not being shown here). Finally, the 
labelled dotted lines are contours of constant age as indicated. The shaded region shows the parameter space not excluded at greater 
than 95 per cent on any single piece of data. 



are taking the dark matter to be cold, it is a reasonable 
hypothesis that the total density of these systems is bigger 
by a factor flo/fln, where J7b = 0.0125ft -2 is the average 
baryon density given by nucleosynthesis. If M is the typical 
mass of the systems, this implies that the fraction /(> M, z) 
of the total mass that resides in bound objects of mass at 
least M at redshift z = 4.0 satisfies 



Therefore we take the minimum m ass of damped Lyman 



/(> M, z = 4.0) > (0.088 ± 0.024) ft 



Qqz 



1 + z 



(32) 



where a 20 per cent uncertainty in the baryon fraction has 
been added in quadrature to the observational uncertainty. 

Since we want a lower bound on the density perturba- 
tion we take the lower end of the error bar. Bearing in mind 
that there is no evidence that damped Lyman alpha systems 
at high redshifts are completely collapsed objects, as we only 
observe their baryonic component which is able to collapse 
faster through radiative cooling (e.g. Katz et al. 1994), we 
conservatively assume that these systems are more akin to 
collapsing protospheroids (see also Lanzetta et al. 1995). In 
order to reflect this choice we will use 8 C = 1.5 in the Press- 
Schechter calculation, which some numerical studies (e.g. 
Monaco 1995) have shown is associated with the time-scale 
of gravitational collapse of a perturbation along its first two 
collapsing axes, i.e. 'filament' formation. Also, in order to be 
compatible with lower redshift observations, the collapsing 
protospheroids have to be massive en ough eventually to giv e 
rise to rotationally supported discs ( Lanzetta et al. 1995 ). 



alpha systems to be 10 u ft M Q (Haehnelt 1995), which 



corresponds to a circular velocity of about 75 kms -1 . Al- 
though formally the constraint as calculated above is only a 
la lower limit, it is almost unchanged by going to 2a. 



5 DISCUSSION 

We plot the data that we have discussed in two separate 
ways. The first is direct contouring of the observations in 
the fio-/i plane, in Fig. 1. For this figure we have fixed the 
normalization of the spectrum by the COBE measurement, 
taking advantage of its small error bar. It turns out that 
the constraint based on the abundance of damped Lyman 
alpha systems is very weak (both in the critical density case 
and for general Qo) as compared to other constraints, and 
so for clarity we do not plot it. All other data play some role 
in constraining the allowed parameters, though the shape 
parameter and cluster abundance allow very similar regions. 
We use the age constraint to cut off the region at the very 
conservative value of 10 Gyr. 

The second type of plot, Fig. 2, shows chi-squared con- 
tours of the data. The only difference in input data is that we 
treat the COBE data with their uncertainty, so that at each 
point in parameter space the chi-squared statistic is that for 
the optimal normalization. The chi-squared plot has the ad- 
vantage of producing a simple summary of the constraints, 
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Figure 2. A chi-squared plot in the Qo—h plane. This plot includes all the data that we have discussed. In contrast to Fig. 1, the COBE 
data are included with their uncertainty error bar. The main compromise is the condensation of the Peacock & Dodds data set down to 
just the shape parameter; otherwise its slightly low chi-squared is excessively generous to the other data as discussed. The inner contour 
corresponds to 68 per cent confidence, the outer one to 95 per cent confidence. Age contours are also shown, but are not included in the 
chi-squared. The chi-squared plot supports the picture formed from Fig. 1, with the drawback of concealing which data are primarily 
responsible for the trends. 



but the drawback that one cannot tell which of the data are 
predominant in contributing to the constraints. 

Most of the recent literature on structure formation 
models has concentrated attention either on retaining firj = 
1 and making other modifications such as introducing a hot 
dark matter component, or on reducing Qo all the way down 
to 0.2 or 0.3. We notice that the best fits with the new COBE 
normalization favour rather higher values, the lowest permit- 
ted being about Qo — 0.35. Further, good fits are available 
for the whole continuum of Qo values above this, for a suit- 
able choice of Hubble parameter. Without inputting extra 
information on the preferred values of h, the observational 
data indicate no particular preference for any value of Qo- 

A variety of recent measurements of the Hubble param- 
eter have favoured higher values (e.g. Schmidt et al. 1994; 
Freedman et al. 1994). While we feel that the situation has 
yet to be completely closed, it is interesting to examine the 
reduction in parameter space implied by choosing h > 0.6. 
This still allows a fit to all the data (even allowing an age 
over 12 Gyr), but such a constraint requires that Qo falls in 
a narrow band between 0.35 and 0.55. 

So far in this paper we have assumed that the spectral 
index of the primordial curvature perturbation is n — 1. In- 
flation models typically predict some degree of 'tilt' in the 
spectrum, so that n is not precisely 1. The degree of tilt pre- 
dicted by inflation is highly model-dependent, ranging from 
negligib le up to a few tenths depending on the inflationary 
model ( Liddle fc Lyth 1993 ). There is some preference for 



tilting to n < 1 but models also exist that act in the op- 
posite direction. There is some reason to believe that tilt 
is rather more likely in open inflationary models, since spe- 
cial physics is being invoked on scales around the curvature 
scale. We end with a short discussion of the effect of tilt. 

In the case of critical density, there is a strong desire 
to remove short-scale power from the spectrum, which both 
tilt and gravitational waves are capable of doing. For low 
densities, the spectrum has already had its shape altered 
by the shifting of matter-radiation equality, and so there is 
less freedom to shift n from the scale-invariant value. In the 
absence of a definite prediction of tilt from an inflationary 
model, we shall examine the two cases n — 0.9 and n = 1.1. 
The effect of gravitational waves in an open universe has not 
been successfully quantified yet and we do not include them 
here. 



Tilting corresponds 
dependent, given by 



to taking Sjj(Qo) to be scale- 



5 2 H {k, Qo) = S 



(33) 



where kio is some normalization scale where the normalized 
spectra at different n are assumed to cross (for a given Qo)- A 
detailed normaliz ation of tilted open models along the lines 
of Gorski et al. ( 1995 ) has not been provided so we need 
some improvization to normalize the models. We assume 
that the tenth microwave anisotropy multipole, which acts 
as a pivot point for the COBE data, is unchanged by tilt; 
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hence the notation kio which is the effective scale of the 
tenth multipole. The prefactor 5jj(Oo) is the normalization 
for n = 1, given by equation We find the scale fcio from 
the 'window function' which describes how different scales 
contribute to the tenth multipole; we take fcio to be the 
scale at which the window function, calculated usin g onl y 
the Sachs-Wolfe effect as in Garcfa-Bellido et al. (1995), 
peaks. We find that fcio is very well fitted by the surprisingly 
simple relation fcio = 6floaoHo- With this, we can then use 
equation (^) in equations ([j]) and (Q). This improvization 
should work very well until £Iq becomes smaller than about 
0.3. 

All other data remain the same, though we now need a 
more general expression f or the shape parameter , which at 
95 per cent confidence is (Peacock & Dodds 1994) 



T = 0.255Z 



+ 0.32 



(34) 



The results are shown in Fig. 3. Increasing n has the 
effect of shifting the allowed band to lower fio , and decreas- 
ing n shifts it to higher Q , roughly in accordance with 
Af2o — —An/2. Assuming the range 0.9 < n < 1.1, we 
therefore find that the width of the band in the £lo~h plane 
is increased, so that for h > 0.6 the allowed range is roughly 
0.30 < Q < 0.60. 

In conclusion, we have made a thorough investigation 
of linear theory constraints on cold dark matter models 
in genuinely open universes, on the assumption that the 
spectrum of the primordial curvature perturbation is scale- 
independent. We have also placed these models in their in- 
flationary cosmology context. The normalization to COBE 
provided by GRSB allows a much more precise comparison 
with observations than has been made previously. We have 
included a treatment of the abundances of both clusters and 
damped Lyman alpha systems; although these have proved 
constraining for various types of model such as mixed dark 
matter models, they are easy to satisfy here. On the whole, 
the new constraints that we have computed support the 
allowed parameter space from earlier considerations rather 
than reduce it. We conclude that there is a substantial pa- 
rameter space still viable for these models. 
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